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Abstract 

Background:The mechanisms behind heterologous immunity and non-specific effects of vaccines on mortality are not well 
understood. We examined associations between cytokine responses and subsequent mortality in low-birth-weight infants in 
Guinea-Bissau. 

Methods: A low-birth-weight trial randomized children to Bacille Calmette-Guerin (BCG) at birth or later according to local 
policy. Blood samples were obtained from a sub-group at age 6 weeks. Interleukin (IL)-5, IL-10, IL-13, interferon (IFN)-y, and 
tumor necrosis factor (TNF)-a were measured in whole-blood cell cultures stimulated with lipopolysaccharide (LPS), 
phytohaemagglutinin (PHA), or purified protein derivative (PPD). The outcome was mortality between bleeding and 1 year 
of age. Non-linear associations between cytokine responses and mortality were examined. 

Results: Cytokine measurements were available from 390 children. The mortality rate (MR) was high (6.8/100 person-years- 
observation (PYO)). Both low and high cytokine responses to LPS and PHA were associated with high mortality (MR up to 
25/100 PYO in the lowest 10% and 9.2/100 PYO in the highest 10%). In BCG-vaccinated children, higher IFN-y responses to 
PPD were associated with better survival (MR ratio = 0.43 (0.24-0.77)). 

Conclusions: Data presented a rare opportunity to explore associations between cytokine responses and mortality. Both 
low and high cytokine responses were associated with high mortality; a balanced response to invading pathogens seems 
preferable. 
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Introduction 

Although infant mortality has declined in low-income countries 
over recent decades, the level is still unacceptably high. A large 
part of the decrease in mortality can be ascribed to vaccinations. 
Recent research has found that vaccines - in addition to the 
specific disease-preventive effects - also have unexpected general 
effects [1-5]. These so-called non-specific or heterologous effects 
are often beneficial by decreasing mortality from other infectious 
diseases than the target disease. Little is known about the 
biological mechanisms behind the non-specific effects, but since 
they are observed in areas where mortality is caused by infectious 
diseases, changes on the immune system may be involved. 

Development of user-friendly immunological methods has made 
it possible to study the immunological effects of vaccines in 



population-based studies. At the Bandim Health Project in 
Guinea-Bissau, we have conducted several randomized trials to 
assess the impact of vaccines on overall mortality, and in 
subgroups of trial participants we have studied the effect of 
vaccines on immunological markers [6-10]. To understand the 
mechanisms behind the non-specific effects of vaccines on overall 
mortality, it is likewise important to identify immunological 
markers of subsequent mortality. Data to explore such associations 
are rare. In the present study we linked results of blood sample 
measurements with data on subsequent mortality to identify 
potential immunological markers of mortality. The association was 
studied in a cohort of low-birth-weight (LBW) children with high 
mortality. 
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Materials and Methods 

Setting 

The Bandim Health Project (www.bandim.org) has implement- 
ed a health and demographic surveillance system (HDSS) which 
covers a population of approximately 100,000 people in six 
suburban districts of Bissau, the capital of Guinea-Bissau. 

Participants 

We used data from the randomized trial "Should Low Birth 
Weight Infants Be Vaccinated With BCG Vaccine at Birth in 
Developing Countries?" (ClinicalTrials.gov, NCT00146302) [5]. 
Bacille Calmette-Guerin (BCG) vaccination was given at birth to 
neonates with a low-birth-weight (LBW) below 2.5 kg. From 2002, 
LBW children born at the national hospital or at two health 
centers in the study area were randomized to receive early BCG (+ 
BCG) or to follow the local policy to receive BCG when the child 
had gained weight, often together with DTP and OPV scheduled 
at 6 weeks of age (—BCG). All children were followed up with 
home visits at 2, 6, and 12 months of age. Causes of death were 
assessed by verbal autopsy. 

Processing of Blood Samples and in vitro Stimulations 

From June 2004— July 2005 we collected blood samples from a 
subgroup of participants when they had reached 6 weeks of age. 
The blood samples, approximately 600 uL, were taken by finger 
prick into heparinized tubes and kept at a maximum of 4 hours at 
ambient temperatures until processing. The whole-blood assay 
and cytokine measurements have been described in detail 
previously [6]. In brief, heparinized whole blood was diluted 
1:10 with RPMI-1640 (Invitrogen, Groningen, Netherlands) 
supplemented with 2 mmol glutamate/L, 1 mmol pyruvate/L, 
100 IU penicillin, and 100 |xg streptomycin/mL. Cultures were 
established in a total volume of 200 uL in 96-well U plate (Nunc, 
Roskilde, Denmark). Supernatants were collected in cultures 
stimulated with lipopolysaccharide (LPS) (100 ng/ml; Sigma- 
Aldrich chemie, Zwijndrecht, the Netherlands) (1 day incubation), 
stimulated with purified protein derivative (PPD) of Mycobacterium 
tuberculosis (10 ug/mL; Statens Serum Institut, Denmark), phyto- 
haemagglutinin (PHA) (2 ug/ ml; Wellcome Diagnostics, Dartford, 
UK) and tetanus toxoid (TT) (1.5 Lf/ml; NVI, Bilthoven, the 
Netherlands) (3 days incubation) and in medium-only-control 
cultures incubated 1 day (CI) and 3 days (C3). The samples were 
stored at — 40°C in Bissau and transported to Leiden, The 
Netherlands on dry ice. All the samples were analyzed in 
November 2005. Concentrations of interleukin (IL)-5, IL-10, IL- 
13, interferon (IFN)-y, and tumor necrosis factor (TNF)-oc were 
measured simultaneously by the commercial Luminex cytokine kit 
(Luminex Corporation, Austin, TX) and buffer reagent kit 
(BioSource, Camarillo, CA) and run on a Luminex- 100 cytometer 
(Luminex Corporation), equipped with StarStation software 
(Applied Cytometry Systems, Dinnington, United Kingdom). 
The lower detection limits (DL) of the assays were 3, 5, 10, 5 
and lOpg/ml for IL-5, IL-10, IL-13, IFN-y and TNF-a, 
respectively. Concentrations below the DL were denoted as non- 
detectable (ND). The samples were measured in two rounds using 
two different batches. 

Statistical Analysis 

All analyses were performed using Stata 12 (Statacorp LP, 
College Station, TX). A significance level of 0.05 was used for all 
analyses. 

For TT stimulations, the cytokine concentrations were very 
similar in stimulated and control samples indicating little effect of 



the stimulation. For all cytokines from the C3 samples (except 
TNF-a) more than 50% of the measurements were ND leaving 
very little information to assess the association with mortality. 
Therefore, these results were not presented. One child who died 
had a missing measurement of IL-13 to PPD and hence the 
mortality analysis using this variable as covariate did not include 
that child. 

The association between cytokine production and mortality was 
analyzed by standardized z-score variables. The cytokine log2- 
concentrations were standardized by BCG randomization, sex and 
measurement batch (subtracting the mean and dividing by the 
standard deviation) creating log-z-score variables with mean 0 and 
variance 1. 

Mortality was assessed at 1 year of age. The follow-up time was 
censored if the family moved and the child was lost to follow-up. If 
a child died, a medical doctor conducted a verbal autopsy with the 
mother or relatives. Crude mortality rates (MR) were presented 
according to percentiles of the log-z-scores. To examine whether 
very low and very high concentrations were associated with high 
mortality, the distributions were split at the percentiles: 1-10%, 
11-50%, 51-90% and 91-100%. 

The association between the mortality rate and a log-z-score 
was analyzed in a Cox-model with time since bleeding as 
underlying time and with age and the log-z-score as covariates. 
ND cytokine concentrations were only known to be below the DL 
(censored at DL); the true values were not known. Hence, the 
corresponding covariate-values in the Cox-regression were also 
censored. Substituting the unknown quantities with a single value 
(e.g. DL) would likely distort the regression. Instead, we used 
multiple imputations (MI) to handle this censoring problem (see 
further details in Methods SI). V-shaped associations between the 
log-z-score and the log-MR were estimated as two log-linear 
relationships below and above the median, restricting the two lines 
to connect at the median (Methods SI). More flexible shapes were 
illustrated graphically. 

We applied complicated methods of analysis to a dataset with 
only 22 deaths. The results were therefore sensitive to model 
assumptions. To accommodate these concerns and to give more 
simple and transparent illustrations of the associations, we also 
illustrated the raw cytokine log 2 -concentrations for dead and alive 
children in Figures SI, S2, and S3, we presented more crude 
analyses stratifying the log-z-scores by percentiles (Table 1 and 2) 
and basic information about each of the dead children was 
presented (Figure 1). The proportional hazards assumption of a 
constant effect across time was explored in an analysis stratified by 
time since bleeding (Table 3). 

Results 

Cytokine concentrations were available from 390 children in the 
LBW trial. The majority of the children (88%) were enrolled in the 
main trial before 7 days of age. The median age at blood sampling 
was 40 days (IQR: 36-45). A total of 199 (51%) children received 
BCG at birth, 1 7 of these received DTP before blood sampling. 
Among the 191 children randomized to not receive BCG at birth, 
24 received BCG only before blood sampling, 9 received DTP 
only, and 14 received both BCG and DTP. OPV was always given 
with DTP. During the follow-up to 1 year of age, 22 children died, 
resulting in a mortality rate (MR) of 6.8/100 person-years of 
observation (PYO). The causes of death determined by verbal 
autopsy were listed in Figure 1. Of the 22 deaths, 13 were among 
female children and 1 1 had been randomized to BCG at birth. 
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Figure 1 . Information on children who died: diagnosis, date of death, time from bleeding to death in days, sex, randomization to +/ 
— BCG, vaccination with BCG, DTP and measles (MV), and cytokine log-z-score percentiles. Note: Rl: Respiratory infection; Heart: Heart 
condition; TB: tuberculosis. The percentile indicates which 10% percentile the corresponding log-z-score belongs to. Vaccine information may be 
incomplete since the vaccination card was often not seen after death and children could have been vaccinated between the last home visit and the 
time of death. 

doi:1 0.1 371 /journal.pone.0093562.g001 



Association between Cytokine Responses and Mortality 

Distributions of the raw cytokine log2-concentrations for dead 
and alive children were presented in Figure SI, S2, and S3. As 
seen in Figure 2, the shape of the association between the log-z- 
score of cytokine responses to LPS and PHA on one hand, and the 
mortality rate on the other hand, were generally well approxi- 
mated by V-shaped relationships, except for IL-10 to PHA. 

Control Responses 

Regarding the TNF-ot concentration in the control sample at 
day 1 (CI), the mortality rate ratio (MRR) was 0.81 (95% 
Confidence Interval: 0.35-1.86) below the median and 1.94 (0.98- 
3.83) above the median (Table 1). Hence, a unit increase in the 
log-z-score was associated with 19% lower mortality below the 
median while it was associated with 94% (borderline significant) 
higher mortality above the median. 

LPS Responses 

For the TNF-ot response to LPS, the MRR was 0.62 (0.43-0.90) 
below the median and 1.64 (0.60-4.48) above the median 
(Table 1). The lowest concentrations (i.e. log-z-scores) were 
associated with higher mortality compared with the highest 
concentrations (Figure 2). Among the 10% lowest concentrations 
the mortality rate was 13/100 PYO compared with 6.1/100 PYO 
among the 10% highest concentrations (Table 1). There were 6 
children with a very low TNF-ot response to LPS (Figure 3). Two 
of these 6 children (33%) died. For the IL-10 responses to LPS, the 
MRR was 0.70 (0.44-1.1 1) below the median and 1.68 (0.56-5.07) 
above the median. Thirty -five children (9%) produced no or very 
little IL-10 (Figure 3). Four of these children (11%) died. 



PHA Responses 

The association between TNF-ot to PHA and subsequent 
mortality was similar to that seen for TNF-ot to LPS; the lowest 
concentrations were associated with higher mortality compared to 
the highest concentrations. Among the 10% lowest concentrations 
the MR was 13/100 PYO. For IFN-y, IL-5 and IL-13, the log- 
linear relationships were significantly different below and above 
the median (Table 1). Among the 10% lowest concentrations the 
MRs were higher than 20/100 PYO. The highest concentrations 
were also associated with high MRs (Table 1). 

PPD Responses 

All cytokine responses to PPD were significantly higher among 
the children randomized to receive BCG at birth (+BCG) 
compared to those randomized to receive BCG later (—BCG). 
The associations between mortality and the cytokine responses to 
PPD were therefore analyzed separately for the two groups. In the 
— BCG group, there were more than 50% ND concentrations for 
IFN-y, IL-5 and IL-13 and there was no association between 
TNF-ot and IL-10 responses and mortality (data not shown). In the 
+BCG group the log-linear relationships were the same below and 
above the median (Figure 2). There was a significant decreasing 
association between mortality and IFN-y; modeled as a single log- 
linear association across the whole range of concentrations (not 
split at the median) the MRR was 0.43 (0.24-0.77) (Table 2). The 
association also tended to be decreasing for IL-5 (MRR = 0.55 
(0.29-1.05)) and IL-13 (MRR = 0.71 (0.38-1.30)). 

BCG Vaccination 

Though not significant, the concentrations of TNF-ot and IL-10 
were larger in the +BC G group compared with the — BCG group for 
C 1 , LPS and PHA. The IFN-y response to PHA was also larger in + 
BCG, while there were no differences for IL-5 and IL- 1 3. The overall 
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Table 2. The associations between the mortality rate (MR) and the log-z-scores of the cytokine responses to PPD were examined 
in children randomized to BCG at birth (+BCG). 





MR/100 PYO (Deaths) 




%ND 


Percentile 








MRR 


1-10 


11-50 


51-90 


91-100 


Overall 


+BCG 


TNF-a 


2% 


6.0 (1) 


7.8 (5) 


3.0 (2) 


20 (3) 


0.99 (0.54-1.78) 


IL-10 


1% 


5.9 (1) 


4.6 (3) 


9.2 (6) 


6.4 (1) 


1.24 (0.78-1.96) 


IFN-y 


6% 


28 (4) 


7.6 (5) 


3.0 (2) 


0.0 (0) 


0.43 (0.24-0.77) 


IL-5 


29% 


21 (5) 


5.3 (3) 


4.5 (3) 


0.0 (0) 


0.55 (0.29-1.05) 


IL-13 


18% 


12 (2) 


6.3 (4) 


6.1 (4) 


0.0 (0) 


0.71 (0.38-1.30) 



Notes: Crude mortality rates (MRs) were presented according to percentiles of the cytokine log-z-scores. The associations between the mortality rate and the log-z- 
scores were the same below and above the median. An overall mortality rate ratio (MRR) was therefore presented. 
doi:1 0.1 371/journal.pone.0093562.t002 



associations between cytokine responses in C 1 , LPS and PHA and 
subsequent mortality were similar in the +BCG group and the —BCG 
group (data not shown). 

Timing of Death 

The higher mortality associated with high TNF-a in CI was 
mainly observed 0-4 months after bleeding (Table 3). The higher 
mortality associated with low responses of TNF-a to LPS or low 
responses of TNF-a, IFN-y, IL-5 and IL-13 to PHA were mainly 
observed 0-4 months after bleeding. On the other hand, high 
responses of TNF-a, IFN-y, IL-5 and IL-13 to PHA were mainly 
associated with high mortality>4 months after bleeding. 

Discussion 

Main Findings 

Mortality among low-birth-weight children in Guinea-Bissau is 
high, in the present cohort the mortality rate was 6.8/100 person- 
years of observation. The immunological study offered a rare 
opportunity to explore associations between cytokine responses 
and subsequent mortality risk. Both low and high cytokine 
responses to stimulation with LPS and PHA were associated with 
high mortality during follow-up. Deaths occurring shortly after 
bleeding were mainly associated with low responses, while deaths 
occurring longer after bleeding were associated with high 
responses. In BCG-vaccinated children a higher IFN-y response 
to PPD was associated with better survival. 

Limitations 

Only children surviving to the age of 6 weeks could be included 
in the immunological study. Thus, the observed associations 
between cytokine responses and subsequent mortality only applied 
to surviving children. We measured stimulated cytokine produc- 
tion from whole blood. Pathogens are often fought in less 
accessible compartments as lymph nodes or mucosa associated 
lymphoid tissue and whole blood may not reflect the capacity in 
these compartments. Infants who died shortly after bleeding may 
have had an ongoing infection. Indeed, we found increased 
mortality risk in children with a high TNF-a concentration in the 
control sample. Peripheral mononuclear blood cells could 
therefore have exited circulation and homed to the secondary 
lymphoid organs. Unfortunately, we do not have blood cell counts 
to corroborate this. Hence, it is a concern that the cause of death 
caused low cytokine production rather than the other way around. 



However, Figure 1 shows that children who had low responses to 
one stimulant generally had high responses to other stimulants. 
Thus, there were some cells present in the blood sample to 
produce cytokines. Overall, the study was small and explorative 
and the associations should be interpreted with appropriate 
caution. 

Interpretation 

Stimulation of whole blood is a measure of the potential for 
activated cytokine production. LPS from gram-negative bacteria 
binds to LPS-binding protein and this complex binds to CD 14, 
which is predominandy found on monocytes or macrophages. 
CD 1 4 signals through Toll like receptor (TLR)-4, leading to NF- 
kB activation and the production of TNF-a and IL-10 [11,12]. 
Whole blood cytokine production after LPS stimulation was found 
to correlate with monocyte production [13] and thus serves as a 
measure of activation of the innate immune response. Conversely, 
PHA is a lectin that stimulates T-cells [14]. The subsequent 
response may be interpreted as the potential for cytokine 
production by cells partaking in the adaptive immune response. 
IFN-y is critical for innate and adaptive immunity against viral 
and intracellular bacterial infections [15]. IL-5 and IL-13 are 
classical Th2 cytokines, promoting B cell growth and maturation, 
activating eosinophil granulocytes, promoting the defense against 
parasites, in addition to playing active roles in asthma and allergy 
disorders [16]. 

In general, mortality was associated with low but also with high 
cytokine responses to LPS or PHA. We hypothesize that the high 
mortality among children with a low potential for cytokine 
production after stimulation of cells from both the innate and 
adaptive immune system reflects a weakened potential to fight 
intruding pathogens. The condition could be compared to 
immune paralysis, a term often used to describe the increased 
susceptibility to secondary infections after sepsis [17] or the 
immune exhaustion among patients with untreated HIV [18-20]. 
It was recently reported that an indeterminate QuantiFERON-TB 
test, due to a very low positive control in terms of the IFN-y 
response to PHA, was associated with high mortality within the 
subsequent 6 months of follow-up [21]. This supports that a low 
IFN-y response to PHA is a predictor of subsequent mortality. It is 
interesting that very high cytokine responses also were associated 
with increased mortality. We speculate that this reflects a higher 
risk of an exaggerated, harmful immune response among these 
individuals. Further studies are required to explore this. 
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Figure 2. Flexible shapes of the associations between the cytokine log-z-scores and the mortality rate. Notes: The curves show mortality 
rate ratios comparing the mortality rate of children with a given log-z-score and age to the mortality rate of children aged 38 days at bleeding with a 
log-z-score equal to zero (Methods SI). For PPD the association is only shown in the group randomized to BCG at birth (+BCG) as there was little 
cytokine production induced by PPD in the group randomized to BCG later (—BCG). 
doi: 10.1 371 /journal.pone.0093562.g002 



Cytokine Production in Response to PPD 

Among the children randomized to receive BCG at birth, there 
was a significant negative association between the IFN-y response 
to PPD and subsequent mortality. The findings corroborate the 
observation that among BCG vaccinated children, developing a 
BCG scar and a large PPD response in vivo is associated with 



reduced all-cause mortality [22]. It remains to be tested whether a 
strong response reflects an underlying stronger immune system or 
whether the association reflects that the children who were 
vaccinated correctly [23] had better survival, i.e. a beneficial non- 
specific effect of BCG. 



LPS, TNF-a 



LPS, IL-10 
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Figure 3. Histograms of the log 2 -concentrations of TNF-a and IL-1 0 in response to stimulation with LPS. Notes: The distribution of TNF- 
cx showed a tail of low concentrations and 6 outlying observations. The distribution of IL-10 included a mixture of responders with normally 
distributed measurements and non-responders (12%) with no or a very weak response to LPS. 
doi: 10.1 371 /journal.pone.0093562.g003 
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Perspectives 

Our results support that a well-functioning state of the immune 
system depends on a subtle balance between strength and 
moderation in responsiveness to foreign pathogens. Specifically, 
we showed that both the upper and lower extreme cytokine levels 
were associated with an increased mortality risk. Several studies 
have explored associations between cytokines and sudden infant 
death [24-28]. One study concluded that very high cytokine 
concentrations in serum or cerebrospinal fluid, obtained during 
autopsy, were a possible explanation for the cause of death [27]. A 
study of extremely low-birth-weight children found bronchopul- 
monary dysplasia or death to be associated with higher 
concentrations of IL-ip, -6, -8, -10 and IFN-y and lower 
concentrations of IL-17 [24]. In this study the cytokine concen- 
trations where analyzed as linear risk-factors; the risk was not 
assessed specifically in the lower and upper extremes. The 
association between anergy and morbidity has been investigated 
in a few studies. Anergy is a state of non-responsiveness to 
stimulation in for instance skin-tests. It may be caused by low 
inflammatory or high regulatory cytokine responses locally [29]. A 
study from Bangladesh reported a higher frequency of diarrheal 
morbidity in children with lack of a delayed type hypersensitivity 
(DTH) response to recall antigens applied by skin test [30]. 
Similarly, a strong associationn was found between the lack of a 
DTH response and follow-up mortality in hospitalized Canadian 
adults [31]. However, the studies did not analyze the association 
specifically in the upper extreme strata. The present study 
indicates that in vitro cytokine responses may serve as useful 
surrogate markers of health status, at least in low-birth-weight 
children, and hence may have potential to serve as outcomes in the 
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